arXiv:1504.05585vl [astro-ph.GA] 21 Apr 2015 


Accepted to ApJ, March 14, 2015 

Preprint typeset using XAT^X style emulateapj v. 5/2/11 


AGES OF STAR CLUSTERS IN THE TIDAL TAILS OF MERGING GALAXIES 

A. J. Mulia^, R. Chandar^ B. C. Whitmore^ 

Accepted to ApJ, March I 4 , 2015 

ABSTRACT 

We study the stellar content in the tidal tails of three nearby merging galaxies, NGC 520, NGC 
2623, and NGG 3256, using BVI imaging taken with the Advanced Gamera for Surveys on board the 
Hubble Space Telescope. The tidal tails in all three systems contain compact and fairly massive young 
star clusters, embedded in a sea of diffuse, unresolved stellar light. We compare the measured colors 
and luminosities with predictions from population synthesis models to estimate cluster ages and find 
that clusters began forming in tidal tails during or shortly after the formation of the tails themselves. 
We find a lack of very young clusters (< 10 Myr old), implying that eventually star formation shuts 
off in the tails as the gas is used up or dispersed. There are a few clusters in each tail with estimated 
ages that are older than the modeled tails themselves, suggesting that these may have been stripped 
out from the original galaxy disks. The luminosity function of the tail clusters can be described by 
a single power-law, dNjdL oc with —2.6 < a < —2.0. We find a stellar age gradient across some 
of the tidal tails, which we interpret as a superposition of 1) newly formed stars and clusters along 
the dense center of the tail and 2) a sea of broadly distributed, older stellar material ejected from the 
progenitor galaxies. 

Subject headings: galaxies: star clusters: general, galaxies: interactions, galaxies: star formation, 
galaxies: star clusters: individual (NGC 520), galaxies: star clusters: individual 
(NGC 2623), galaxies: star clusters: individual (NGC 3256) 


1. INTRODUCTION 

Tidal tails form during the interaction between galax¬ 
ies, and are composed of stars, gas, and dust ejected from 
the original galactic disks. These interactions are often 
accompanied by one or more bursts of star and cluster 
formation in the main bodies of the parent galaxies, ac¬ 
cording to simulations (e.g., Mihos, Bothun, & Richstone 
1993) and observations of star clusters (e.g., Whitmore 
& Schweizer 1995). 

The star cluster formation processes that occur within 
tidal tails, however, are not well understood. For ex¬ 
ample, building on pioneering work by Knierman et al. 
(2003), Mullan et al. (2011; hereafter Mil) used V and 

1 band photometry from Wide Field Planetary Camera 

2 (WFPC2) on the Hubble Space Telescope (HST) to 
measure the brightness and colors of point-like sources 
both in and out of the tail region for 23 tidal tails. They 
found a statistically significant number of clusters in 10 
out of their 23 tails, and they concluded that the pres¬ 
ence of tail clusters depends on several factors, including 
tail age, H I density, and surface brightness of the tail. In 
this work we revisit two of these tails using higher qual¬ 
ity data taken with the Advanced Camera for Surveys 
(AGSjSon HST/Wide Field Camera (WFC) and detect 
clusters in both of them for the first time. 

Even less is known about the ages of star clusters in 
tidal tails, de Grijs et al. (2003) found ^ 40 clusters each 
in the Tadpole galaxy with rough ages of 175 ± 25 Myr 
and in the Mice galaxies (NGC 4676) with ages ^ 100T20 
Myr. Most of these clusters appear to have formed in the 
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tails since their estimated ages are younger than those of 
the tails; de Grijs et al. estimated a dynamical age of 
400 — 800 Myr for the Tadpole, and Privon et al. (2013) 
suggested a formation time of ^ 175 Myr ago for the 
Mice. Tran et al. also observed the Tadpole galaxy, and 
found 42 extremely young tail clusters (^4 — 5 Myr; 
2003). Bastian et al. (2005) found young massive star 
clusters in the tails of NGC 6872 and was able to esti¬ 
mate ages and masses for individual clusters using UBVI 
and Ha photometry. No star clusters older than the dy¬ 
namical age of the tidal tail were found (at least down 
to their completeness limit of My = —10.4), once again 
suggesting that the clusters formed from the tidal debris. 
They also found a large population of clusters that are 
very young (< 10 Myr). While looking for spatial trends, 
Bastian et al. found that older tail clusters tend to be 
located closer to the main bodies, while younger clusters 
are spread throughout the tail. They argue that these 
findings suggest an initial burst of cluster formation as 
the tails form. After this burst, as the tails expand, the 
gas clouds cool and condense, eventually becoming sites 
of individual cluster formation. 

While tidal tails provide a unique and interesting 
environment to study star clusters, they also provide 
key information for improving simulations of interact¬ 
ing/merging galaxies. Efforts to reproduce mergers using 
simulations focus on tidal tail morphology to constrain 
details of the galaxy interaction. Ages of star clusters in 
tails provide constraints on the merger timescale. The 
spatial distributions of clusters in the main bodies and 
tails helps to discriminate between different prescriptions 
for star formation within the simulations. Eor exam¬ 
ple, dynamical models of the interacting system NGC 
7252 show that density-dependent star formation results 
in a more centrally concentrated distribution of main- 


2 


body clusters, while shock-induced star formation tends 
to spread cluster formation throughout the galaxy, in¬ 
cluding in the tails (Chien & Barnes 2010). 

The goal of this work is to better understand how clus¬ 
ters form and are distributed during the merging of two 
galaxies. We focus our study on three merging systems, 
NGC 520, NGC 2623, and NGC 3256. This paper is 
set up in the following manner. In Section 2 we discuss 
the observations and photometry, as well as cluster selec¬ 
tion. Section 3 presents the cluster densities and cluster 
luminosity functions, as well as colors and ages of clus¬ 
ters and of the diffuse light. In Section 4 we discuss the 
galaxy merger history and interpret our results. Section 
5 summarizes the work and presents conclusions. 

2. OBSERVATIONS, DATA REDUCTION, AND CLUSTER 
SELECTION 

2.1. Observations and Data Reduction 

Observations of NGC 520 and NGC 2623 were taken 
with the WFC on the ACS on HST as part of pro¬ 
gram GO-9735 (PI: Whitmore). NGC 520 was observed 
in 2004 October using the filters F435W (« B in the 
Johnson-Cousins system; exposed for 2400 s), F555W 
1/; 1000 s), and F 8 I 4 W /; 1200 s). The entire 
galaxy, including the majority of its tidal tails, were cov¬ 
ered in a single WFC frame. Observations of NGC 2623 
were taken in 2004 February using the same filters with 
exposure times of 3730, 1206 and 2460 s, respectively. 
It was also covered in a single WFC frame. Images of 
NGC 3256 come from two separate observing proposals. 
F555W filter observations of NGC 3256 were obtained 
in 2003 November as part of the program GO-9735 (2552 
s). WFC observations using F435W and F 8 I 4 W filters 
were taken in 2005 November as part of program GO- 
10592 (PI: Evans) for 1320 and 760 s, respectively. Both 
fields cover most of NGC 3256’s eastern tail, as well as 
its main body. Figure [T] shows the WFC fields of view 
used in this study. 

The raw data were processed through the standard 
ACS pipeline. The reduced, multidrizzled images were 
taken from the Hubble Legacy Archive and have a scale 
of 0.05" per pixel. Figure [2] shows F 8 I 4 W images of our 
targets, with the tidal tail portions studied in this work 
labeled. We will refer to specific tails within a galaxy by 
abbreviating the location of each tail (i.e. the northern 
tail of NGC 520 will be called NGC 520N). We define 
the tail edges as the locations where the counts drop to 
3cr above the general background level. 

Figure [ 2 ] shows that NGC 520 has two tidal tails, com¬ 
monly referred to as the north and south tails. Both of 
these tails extend beyond the field of view in Figure [2l 
Not seen in Figure [2] is the dwarf galaxy UGC 957, at 
about 6' north-west of NGC 520. NGC 2623 is in the 
middle panel of Figure O It exhibits two nearly symmet¬ 
ric tidal tails, although the northern tail is wider than the 
southern tail. An extremely bright, triangle-shaped re¬ 
gion south of the nucleus of the galaxy, which we call the 
pie wedge, is also highlighted. NGC 3256 is on the right 
in Figure [21 While it has two tails (east and west), only 
part of the eastern tail was observed. The three merg¬ 
ers are next to each other on the Toomre Sequence, in 
the order NGC 520, NGC 2623, and NGC 3256 (Toomre 
1977). General properties of these systems are listed in 
Table m 


2.2. Cluster Selection 

Star cluster candidates were found by running the 
IRAlfl task DAOFIND on the I band images. A few 
thousand objects were detected in each galaxy. We per¬ 
form BVI aperture photometry on all sources in NGC 
520 and NGC 3256 using an aperture size of 3 pixels in 
radius and a background area of radii 8 to 11 pixels using 
the IRAF task PHOT. Because of the greater distance 
to NGC 2623, we use a smaller aperture of 2 pixels (with 
the same background area) in order to minimize contam¬ 
ination from nearby sources. In the pie wedge region, we 
estimate the background level using annuli between 5- 
8 pixels, in order to minimize the impact of crowding. 
Magnitudes were converted to the VEGAMAG system 
using zeropoints calculated from the ACS zeropoint cal- 
culatoiEI. Aperture corrections were performed in NGC 
3256 by selecting the brightest ^ 20 isolated clusters in 
each filter and finding their mean 3-10 pixel magnitude 
difference. We then applied the 10 pixel to infinity aper¬ 
ture corrections taken from the encircled energy catalog 
for ACS/WEC in Table 3 of Sirianni et al. (2005). The 
final corrections were 0.45, 0.42, and 0.54 magnitudes 
for R, F, and /, respectively, and they were applied to 
each cluster. Unfortunately, the background level was 
too high in NGC 520 and NGC 2623 for reliable aper¬ 
ture correction determinations. We find that clusters in 
NGC 520 have a similar range of sizes as those in NGC 
3256 and therefore applied the aperture corrections de¬ 
rived from NGC 3256 to the NGC 520 clusters. Clusters 
in the more distant NGC 2623 are nearly indistinguish¬ 
able from point sources; we assume the 2 pixel to infin¬ 
ity aperture corrections from Table 3 of Sirianni et al. 
(2005). 

We measure EWHM values for all detected sources us¬ 
ing the ISHAPE software (Larsen 1999). ISHAPE mea¬ 
sures EWHMs by convolving the point spread function 
with a King profile, then performing a calculation to 
test the goodness of fit to each individual cluster (King 
1966). ISHAPE iterates through different values for the 
effective radius until a minimum is found. 

We select cluster candidates using a combination of au¬ 
tomated selection criteria (listed in Table |2j) and visual 
inspection of each candidate (see Bastian et al. 2012 for 
further discussion on the “hybrid” method). The selec¬ 
tion criteria in my, EWHM, and magnitude error (ay) 
remove most foreground stars and background galaxies. 
We made a cut at the bright end of my = 20, in order 
to eliminate saturated foreground stars (brighter than 
any cluster expected in these mergers). The faint end 
of our cluster catalog likely suffers from incompleteness. 
We restrict our sample to luminosities brighter than the 
point where the luminosity function (discussed in Sec¬ 
tion 3.2) turns down rather than continuing to increase 
in a power-law fashion. We use cuts in EWHM (size) to 
separate stars and background galaxies from clusters, to 
the highest degree possible. In NGC 3256, we exclude 
sources with EWHM values less than 0.2 pixels, because 
the EWHM of the vast majority of likely field stars was 

^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation. 

^ http://www.stsci.edu/hst/acs/analysis/zeropoints 
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Figure 1. 10’ x 10’ images of our galaxy merger sample from the Digitized Sky Survey. Overlaid in each image is the WFC field of view. 
In these and all the following images, north is up and east is to the left. 



Figure 2. ACS/WFC FSIAW images showing the sample of all three galaxies and their tails. 


Table 1 

Galaxy Properties 


Galaxy 

R.A. 

Dec. 

Distance Modulus,^ 

Av 

E(B - V) 

E(V - I) 

NGG 520 

01 24 35.1 

+03 47 33 

32.46 

0.077 

0.025 

0.035 

NGG 2623 

08 38 24.1 

+25 45 17 

34.50 

0.113 

0.036 

0.051 

NGG 3256 

10 27 51.3 

-43 54 13 

32.79 

0.334 

0.107 

0.151 


^Distance moduli come from distances obtained from recession velocity measurements relative to the local group, adopting Hq = 70 km 
s“^ Mpc“^. 


Table 2 

Selection Cuts 


Galaxy 

my 

FWHM 

cry 

NGG 520 

25.0 

0.1 - 2.0 


NGG 2623 

26.5 

< 2.0 

0.1 

pie wedge 

26.5 

< 2.0 

0.2 

NGG 3256 

26.0 

0.2 - 2.0 



measured to be smaller than that value. We also place 
an upper limit of 2.0 pixels in FWHM to exclude very ex¬ 
tended background galaxies. The lower limit on FWHM 
for clusters in NGC 520 was chosen to be 0.1 pixels be¬ 
cause there is a group of young, compact clusters in the 
northern portion of NGC 520S. At the distance of NGC 
2623, clusters are essentially point sources, so FWHM 
measurements cannot distinguish clusters from bright, 
individual stars. Therefore, in addition to an upper limit 
on FWHM of 2.0, we apply a cutoff in magnitude uncer¬ 
tainty, (7y. In the tails, (jy <0.1 did a good job removing 
obvious non-clusters, but in the pie wedge, we relax our 
(Ty cut to 0.2 because the area is more crowded and has 
a very bright background. 

Finally, we visually inspect each cluster candidate. For 


each galaxy we identify high signal-to-noise ration (S/N) 
clusters, as well as obvious foreground stars, as bench¬ 
marks for fainter objects. We choose clusters based on 
their fuzzy appearance, as well as wider radial profiles 
than those of stars. Figure [3] shows three of the brightest 
clusters from each galaxy’s tidal tails, and Figures |4]-[3 
show locations of clusters in each tail. All possible clus¬ 
ters were identified independently of the selection cuts 
listed above, then later compared to the sample obtained 
from the selection cuts. We find excellent agreement be¬ 
tween the two methods in NGG 520 and NGG 3256, and 
we only include sources that both fit our selection cuts 
and visually appear as clusters. Due to the poor S/N 
in NGG 2623, several sources appeared visually as likely 
clusters but had FWHM and/or ay values outside of the 
selection criteria. Because of this low S/N, we rely more 
heavily on visual inspection than selection cuts and in¬ 
cluded some sources with FWHM or ay values outside 
of the cuts. The total number of clusters found in each 
tidal tail {Ndus) is compiled in Tabled and we compile 
the cluster catalog in Table IH 

Distant, red spheroidal galaxies can have a similar 
color and appearance to compact star clusters. To assess 
the possible contamination from background galaxies, we 
examine the colors and locations of sources that made 
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Table 3 

Properties of star clusters. Ttaii [Myr] is the predicted age of the tidal tail as estimated from simulations. A^clus is the total number of 
clusters found in the region. S^us is the star cluster density per kpc^ where My < —8.5. Tdus [Myr] is the approximate age of the bulk 
cluster formation. Median colors ((B - V), (V - I}) are for clusters younger than the age of the tidal tail, a values (from the relation 

dNjdL oc L^) listed come from the aconst-num method. 


Galaxy 

Tail 

Rail 

Nclus 

^clus 

Rlus 

(B-V> 

(V-I> 

a 

NGG 520 

North 

300 

5 

0.00 

- 100 

0.08 ± 0.05 

0.38 ± 0.05 


NGG 520 

South 

300 

93 

0.49 

260 

0.19 ± 0.09 

0.60 ± 0.17 

-2.32 ± 0.36 

NGG 2623 

North 

220 

75 

0.54 

230 

0.10 ± 0.15 

0.48 ± 0.28 

-2.39 ± 0.34 

NGG 2623 

South 

220 

18 

0.32 


0.12 ± 0.13 

0.43 ± 0.17 


NGG 2623 

pie wedge 

< 220 

76 

8.5 

- 100 

0.05 ± 0.15 

0.43 ± 0.19 

-2.02 ± 0.21 

NGG 3256 

East 

450 

141 

0.10 

> 260 

0.21 ± 0.12 

0.58 ± 0.21 

-2.61 ± 0.27 


Table 4 

Cluster catalog. Apparent magnitudes {m\) and photometric errors are listed, as well as individual cluster ages in log years (r). 
Table 0] is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 


Tail 

ID 

RA 

Dec 

UlB 

CTB 

my 

cry 

mi 

cri 

r 

520N 

1 

1 24 31.77 

+3 47 59.03 

24.68 

0.02 

24.01 

0.03 

22.82 

0.02 

9.11 

520N 

2 

1 24 31.02 

+3 48 3.70 

24.63 

0.03 

24.47 

0.04 

24.10 

0.04 

8.06 

520N 

3 

1 24 32.06 

+3 48 6.88 

25.82 

0.07 

24.70 

0.05 

23.41 

0.03 

9.99 

520N 

4 

1 24 29.58 

+3 48 27.96 

24.31 

0.02 

24.23 

0.03 

23.85 

0.03 

8.06 

520N 

5 

1 24 28.77 

+3 48 35.25 

24.41 

0.02 

24.33 

0.03 

23.88 

0.03 

8.06 

520S 

1 

1 24 36.18 

+3 45 37.42 

24.95 

0.03 

24.56 

0.04 

23.76 

0.03 

8.76 

520S 

2 

1 24 36.27 

+3 45 38.32 

25.14 

0.03 

24.40 

0.04 

23.20 

0.02 

9.30 

520S 

3 

1 24 35.59 

+3 45 41.18 

25.00 

0.03 

24.94 

0.05 

24.34 

0.04 

7.59 

520S 

4 

1 24 36.11 

+3 45 43.04 

23.66 

0.01 

23.25 

0.02 

22.46 

0.01 

8.81 



Figure 3. Some of the brightest star clusters in the tidal tails of 
each galaxy. Each image is 25 x 25 pixels. Galaxies are separated 
by column from left to right, in the order NGC 520, NGC 2623 
(tail region), NGC 2623 (pie wedge region), and NGC 3256. 

it through our automated cuts, but are located outside 
of the tidal tail regions, with the reasonable expecta¬ 
tion that most of these are likely background galaxies. 
However, we find no such objects in any of our fields, 
only extended, obvious background galaxies. In addi¬ 
tion, many background galaxies should appear elliptical 
rather than circular. Using IS HA PE, we measure the el- 
lipticity of each cluster candidate and find few elongated 
red sources. We conclude that the fraction of background 
galaxies that have made it through our cluster selection 
criteria is negligible. 

3. RESULTS 

We find a population of star clusters in each tidal tail in 
our sample of galaxy mergers. In this section we present 


the spatial densities, luminosities, and colors of clusters 
in each tidal tail and estimate their ages. 

3.1. Cluster Density 

One quantitative measure of a cluster population is its 
surface number density, which is often used when clus¬ 
ter selection must be done statistically. We measure the 
number density of clusters, Sdus (clusters kpc“^), for 
each tidal tail based on the catalogs in this work. We 
choose a magnitude cutoff of My < — 8.5 (see Section 3.2 
for details on absolute magnitude calculation) to directly 
compare our results to those of Mil (although altering 
the cutoff to My = —8.0 and —9.0 yields qualitatively 
similar results). NGC 520S yielded Scius = 0.49, while 
NGC 520N yielded Hcius = 0.0 because its five clusters 
were all fainter than My = —8.5. The pie wedge in NGC 
2623 possesses the highest cluster density in our sample 
by over an order of magnitude at Sdus = 8.5, while the 
north and south tails yielded Sdus = 0.54 and 0.32, re¬ 
spectively. Lastly, NGC 3256E possessed a relatively low 
density of Edus = 0.10. We find that while all tidal tails 
in our sample contain clusters, there is a large variation 
in Edus for each tail. 

For both NGC 520 and NGC 3256, we are able to di¬ 
rectly compare our values of Edus to those values found 
in Mil, although we note that we cover different areas of 
the tails for both systems. We also note that we adopt 
slightly different distance moduli than Mil, resulting in 
small differences in Edus- Mil found Edus = —0.011 ± 
0.013 in NGC 520N and Edus = 0.004 ± 0.082 for NGC 
3256E, both of which are consistent with zero. While 
our results for NGC 520N agree with Mil, we find a 
clear population of clusters in NGC 3256E brighter than 
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Figure 4. NGC 520 cluster candidates are marked with green circles. On the top is the northern tail. The lower-left panel shows the 
northern portion of the southern tail. On the lower-right is the southern half of the southern tail. 


Mv < —8.5. Mil also measured Hdus foi* My < —7.5 
and My < —6.5, however we do not compare Hcius val¬ 
ues at these cutoffs due to incompleteness. We find that, 
in general, younger tails tend to exhibit higher Sdus^ 
agreement with the results of Mil. 

3.2. Luminosity Functions (LFs) 

The LFs of the tail clusters are shown in Figure [51 
These can be described by a power law, dNjdL oc 
where a is determined from a linear fit to log 
{dN/d{My)). We fit our cluster LFs over the V band 
magnitude range which we believe to be reasonably com¬ 
plete. Magnitudes include distance moduli from Table 
H] as well as aperture corrections, and they have been 
dereddened because of foreground extinction (see Sec¬ 
tion 3.3.1). We are unable to measure the LF in NGC 
520N and NGC 2623S because they contained too few 
clusters. 

We determine the LFs in two different ways, using both 
a fixed magnitude bin width (denoted <aconst-mag) and a 
fixed number of clusters per bin (o^const-num)- The re¬ 
sults from both methods are consistent, and we report 
a values for the latter because it is better at handling 
low number statistics. In all cases, however, we find that 
the two values for a are the same, within the uncertain¬ 
ties. The typical standard deviation among a calculated 
using I^ and V band magnitudes is 0.13, with no 
systematic trend. The LFs in the B and I bands have 
similar slopes to those in the V band. 

In NGC 520S, we find a best fit of a-2.32 ± 0.36 


down to My < —8.0. We restrict our fit to this magni¬ 
tude limit because completeness becomes an issue fainter 
than this. For NGC 2623, we focus on the pie wedge and 
northern tail. We fit the luminosity function in both re¬ 
gions down to My ^ —9.0, and find a = —2.39 ± 0.34 
for the northern tail and a = —2.02 ± 0.21 for the pie 
wedge. The luminosity function in NGC 3256E is best 
fit with a = —2.61 ± 0.27 down to My = —8.0. We find 
no evidence for a deviation from a simple power-law at 
the bright end of the LF in any of our tail regions. 

Several studies have shown that the luminosity func¬ 
tion can range from —1.9 < a < —2.8 for clusters in a 
variety of environments (e.g., Whitmore et al. 2014 and 
references therein). We find that our results for a are 
within this range. 

3.3. Colors and Ages of Clusters in Tidal Tails 

We present the color-color diagrams for clusters in all 
tails in Figure |9l The color coding corresponds to dif¬ 
ferent tails within each merging system. We also show 
with blue symbols the median V — I and B — V colors of 
clusters that apparently formed after the tidal tails (we 
assume tail ages from simulations that are discussed in 
Section 4.1 and compiled in Table El- The solid line in 
each panel in Figure [9] is a stellar population model from 
G. Bruzual & S. Chariot (2006, private communication, 
hereafter BC06; and see also Bruzual & Chariot 2003) 
which predicts the evolution of star clusters from about 
10^ - 10^^ yr. All models assume a metallicity of Z = 
0.02, and a Chabrier (2003) initial mass function. Num- 
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Figure 7. NGC 3256E cluster candidates. 
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NGC 520S 



NGC 3256E 



NGC 2623N 



Pie Wedge 



Figure 8. Luminosity functions of clusters in tails where measurements were possible. Filled circles are derived from the constant-number 
method and are best fit by the dashed line (slope given as aconst—num; see section 3.2 for details). Open circles fit with a dash-dotted line 
(slope given as ^const-mag) are derived from the constant magnitude bin method. 


bers mark the logarithmic age (r) corresponding to the 
population. 

3 . 3 . 1 . The Impact of Internal Reddening on Cluster Colors 

In order to accurately age-date star clusters, we must 
understand how reddening affects their measured colors. 
We correct each cluster for foreground reddening of the 
Milky Way by using Ax values given in the NASA/IPAC 
Extragalactic Database, which are compiled in Table [TJ 
The reddenings were estimated with the standard rela¬ 
tions between them and the extinction. A reddening vec¬ 
tor for Ay = 1 is shown in Figure [9] by the red arrow. 
The reddening vector is nearly parallel to the predicted 
evolution of cluster colors, making it challenging to dis¬ 
entangle the effects of age and reddening in the measured 
colors. For example, a young cluster with some extinc¬ 
tion will have the same colors as an older cluster with no 
reddening. 

Figure [9] shows that the clusters in each tail have a 
measured color spread. In principle, there are a few rea¬ 
sons this spread could exist other than photometric un¬ 
certainties, such as internal (and therefore differential) 
extinction or a range in cluster ages. Below, we examine 
each galaxy individually and determine that internal ex¬ 
tinction does not greatly affect the majority of our cluster 
colors. 

We first focus on NGC 2623, whose color distributions 
are shown in the upper- and lower-right panels of Figure 
[9l The spread in measured colors is the smallest for any 


tail studied here, and instead of spreading along the red¬ 
dening vector, the cluster colors are fairly isotropically 
distributed. This pattern suggests that photometric un¬ 
certainties are primarily responsible for the color spread. 
We confirm this by creating color-color diagrams of NGC 
2623 clusters in different My bins, and while the cluster 
color spread increases with fainter magnitudes, the color 
distribution is centered near the median colors quoted in 
Table O regardless of brightness. If all of the clusters 
in the tails of NGC 2623 formed at the same time, clus¬ 
ters obscured by dust would exhibit redder colors than 
unobscured clusters, causing a color spread along the red¬ 
dening vector. The clusters in both tidal tails of NGC 
2623 appear to have formed over a relatively short pe¬ 
riod of time. This is also true in the pie wedge, where the 
narrow color ranges are also indicative of coeval cluster 
formation. 

NGC 520 and NGC 3256, on the other hand, both ex¬ 
hibit a spread in color along the reddening vector. Here, 
we use the relative locations of blue {B — V < 0.3) versus 
red (H — F > 0.3) clusters, shown in Figures [TOl and [TT| 
to assess whether age or reddening are more likely to be 
responsible for the observed spread in color. In the case 
of NGC 520 (Figure [To]), the bluer clusters tend to trace 
out the central, densest portion of the tail, while the red¬ 
der clusters are more evenly spread throughout the tail. 
Reddening due to dust would be expected to give the 
opposite trend, with reddened clusters mostly near the 
center of the tail. Figure [IT] shows the locations of blue 
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Figure 9. Color-color diagrams. Clusters have been corrected for foreground reddening of the Milky Way by using the Ax values compiled 
in Table [T] Solid black lines are BC06 cluster evolution tracks, and numbers represent log ages of clusters. The reddening vectors are 
shown by the red arrows for Ay = 1- For comparison, the red diamond marks where the tail age (rtaii) would fall on this diagram based on 
the ages shown. Clusters that formed in the tail are younger than the tail and presumably fall above the red diamond on the track. The 
median colors of these young clusters are marked in blue, represented by either a square (northern/eastern tail), triangle (southern tail), 
or a cross (pie wedge). The points in the upper right-hand corners indicate typical errors from photometric uncertainties. 
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Figure 10. Southern tail of NGC 520. Red and blue circles 
indicate clusters that are redder and bluer than our color cut 
{B — V = 0.3), respectively. 

versus red clusters for NGC 3256. No obvious pattern 
is observed, but there is no evidence that red sources 
are primarily found in the brightest, densest portion of 
the tail. These results support the interpretation that 



Figure 11. Eastern tail of NGC 3256. Red and blue circles 
indicate clusters that are redder and bluer than our color cut 
(B — V = 0.3), respectively. 


the larger color spread of the clusters in the tails of NGC 
520 and NGC 3256 is primarily due to a spread of cluster 
ages rather than to varying amounts of cluster reddening. 

We now compare the colors of red tail clusters to clus¬ 
ters in the halos of NGC 520 and NGC 3256, which are 
shown in Figure [121 Selection and photometry of the 
halo clusters are performed using the same methods as 
in Section 2. We conclude that they are fairly similar 
in color, and therefore the red tail clusters are likely to 
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be clusters that formed before the tails but were swept 
out during the tail formation. We will study these halo 
clusters and other main body clusters in more detail in 
a forthcoming paper. All of this evidence supports the 
interpretation that the red sources are older, unobscured 
star clusters. 



Figure 12. Locations of star clusters in the halos of NGC 520 
(left) and NGC 3256 (right). 

3.3.2. Ages 

We obtain ages for each cluster candidate by mapping 
their B—V and V—I colors to the BC06 track and finding 
the closest match (Figure [9]). Ages are fairly accurate in 
the range 10^ < r < 10^ yr. Based on the discussion in 
Section 3.3.1, we assume no internal reddening. 

Figure [9] implies that star clusters in the tails of NGC 
520 and NGC 3256 have a broad range of ages up to 
- 10^-^ yr. The tail ages are marked for each system with 
a red diamond, and clusters with colors below the red 
diamond are older than the tail age and might be old disk 
clusters that predate the merger; we will hereafter refer 
to these as “old” clusters, and we will refer to clusters 
younger than the tail age as “young” clusters. The tail 
clusters in NGC 2623 contain ages up to a few xlO^ yr. 
The lack of observed clusters with ages older than a few 
X10^ yr may simply be because they are too faint to be 
observed in our data at the distance of NGC 2623. 

3.4. Diffuse Light in Tidal Tails 

In addition to distinct stellar clusters, we also study the 
diffuse stellar light in the tails. Because the tidal tails 
form from material stripped from the native disks, this 
diffuse light should consist of a mixture of old disk stars, 
faint newly formed stars, and dissolved clusters. While 
individual stars are unresolved, the integrated colors of 
the diffuse light can help to constrain the age of the stel¬ 
lar populations (e.g., Gallagher et al. 2010). The light 
will be dominated by the brightest and youngest unre¬ 
solved objects, so we expect that our age estimates will 
be lower limits. 

We create color images: B — V^ V — and B — I for all 
three galaxies, and we show the B — I images in Figure 
M The images were boxcar-smoothed with a 4 pixel 
kernel width. The main bodies of the mergers are very 
red because they contain large amounts of dust. The 
tidal tails are easily seen in Figure [13] as bluer regions 
where internal extinction is low (except, perhaps, in NGC 
520N, as also suggested by our analysis of clusters in 
Section 3.3.1). We therefore assume that colors in the 
tail regions are not greatly affected by internal extinction, 
and we believe diffuse light colors trace the ages of stellar 
populations relatively well. 


Table 5 

Diffuse Light Measurements 


Galaxy 

Region 

(V “ I) diffuse 

(B - V) diffuse 

NGC 520 

NA 

1.13 

0.70 


NB 

1.05 

0.60 


SA 

1.14 

0.71 


SB 

0.94 

0.49 


SG 

0.90 

0.47 


SD 

0.88 

0.46 


SE 

0.97 

0.49 

NGG 2623 

NA 

0.96 

0.48 


NB 

1.07 

0.42 


NG 

0.76 

0.33 


ND 

0.90 

0.47 


SA 

0.88 

0.45 


SB 

0.87 

0.40 


SG 

0.92 

0.38 


SD 

0.91 

0.41 


SE 

0.96 

0.42 


P 

0.68 

0.31 

NGC 3256 

EA 

0.77 

0.82 


EB 

0.71 

0.79 


EG 

0.78 

0.39 


ED 

0.83 

0.76 


EE 

0.66 

0.63 


We measure the V — I and B — V colors in different tail 
regions (white boxes in Figure [13] labeled from “A” up 
to “E” for each tail) after subtracting the point sources. 
Colors are calculated from the ratio of fluxes on a pixel 
by pixel basis and corrected for differences in instrumen¬ 
tal zeropoints and foreground extinction. The diffuse 
light regions were chosen to cover most of the optically 
bright portions of the tails while avoiding breaks in field 
coverage and saturated foreground stars. The values in 
Table [S] are the mean color for each region. Figure [14] 
plots these values on a color-color diagram overlaid with 
the BC06 track. Figure [13] reveals that the diffuse stellar 
light in all tidal tails in our sample exhibits a common 
pattern (with the exception of NGC 3256E). While the 
tidal tails are fairly blue {B — I < 1.4), we find that 
there is a gradient with B — I coior; it is redder near 
the edges of the taii and biuer near its axis. As previ- 
ousiy discussed, this is not iikeiy a reddening effect, so 
we believe that we are observing a gradient in the mean 
ages of the stars across these tails. In Section 3.3.1 we 
discussed the tendency for younger clusters to be near 
the center of the tail axis, leaving mostly older clusters 
along the tail edges, especially in NGC 520S. The diffuse 
light images further support a lack of young stars near 
the edges of the tails. We believe this to be the first 
time such a gradient across the tails has been observed 
in these galaxies. 

Eigure [14] reveals that the northern tail of NGC 520 
has the reddest diffuse light of any tail in our sample. 
This is likely due to dust, as the red B — I color in Eig¬ 
ure [13] extends from the main body straight through the 
tail. Region A in the southern tail is also almost cer¬ 
tainly affected by dust. We believe that the colors of 
the diffuse light in the other regions of the southern tail 
are dominated by ages of the underlying populations and 
not extinction because of the absence of dust lanes in this 
tail. These regions have an age r of around 800 Myr. 
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Figure 14. Diffuse light colors for all thre e ga laxies. Regions 
indicated in the legend correspond to Figure fT3l Asterisks mark 
the median color of clusters found in that tidal tail that presumably 
predate the merger. Asterisks inside squares indicate that there 
were very few clusters that predate the merger (< 3). The solid 
black line is a BC06 cluster evolution track. 


r ~ 650 Myr. This would indicate that the diffuse light 
population in region EC is similar to the majority of the 
diffuse populations in the tails of NGC 2623. It would 
be plausible that region EC falls on a local region of en¬ 
hanced star formation, causing its diffuse light colors to 
be bluer than other regions of NGC 3256E. However, re¬ 
gion EC does not include any young star clusters, instead 
containing only one cluster whose colors suggest that it 
predates the merger. The rest of the regions we measured 
are gathered around E — / ~ 0.7 and B — V ~ 0.8 and 
are off the BC06 track. Regardless, these offset colors 
suggest faint population ages of about 800 - 1000 Myr. 
It is likely that most of these diffuse light regions do not 
have colors consistent with the BC06 track because of the 
low surface brightness of NGC 3256E, resulting in large 
photometric errors. It is unlikely that Ha line emission 
could be responsible for the offset from the model, al¬ 
though this emission would move measurements in the 
observed direction relative to the models, because such 
recent star formation is usually quite bright and lumpy. 

4. DISCUSSION 


In NGC 2623, most of the regions in the north tail 
have colors consistent with stellar populations of ^ 700 
Myr old. Region NC, however, is bluer and yields an age 
of ^ 500 Myr. We note that the diffuse light in the pie 
wedge (region P) has B — V and V — I colors close to 
region NC, as seen in Table O which might provide clues 
to its formation, as we discuss further in Section 4.1.2. 

Eigure [TU shows that the diffuse light in NGC 3256 is 
perhaps the most peculiar of our sample. Only region 
EC falls near the BC06 track, with an associated age 


4.1. Comparing Cluster Ages with Tail Ages 

The age of the tidal tails, which we define as rtaii, refers 
to the time since perigalacticon of the two progenitor 
galaxies, because they form from material ejected after 
the galaxies undergo their closest approach. The ages of 
tails can be estimated from simulations and are proba¬ 
bly accurate to within ^ 20% 30 — 100 Myr; Barnes 

& Hibbard 2009). In this section, we briefiy summarize 
previous simulations of our sample of mergers, and we 
compare tail ages from those simulations with our mea- 
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sured cluster ages. 

We find that, for most systems, there is a period of 
increased cluster formation near the time that the tidal 
tails formed. We quantify the time of this cluster for¬ 
mation (hereafter denoted Tdus) by comparing the den¬ 
sity peak of the color distribution with cluster evolution 
model predictions. These density plots are shown in the 
right panels of Figure [151 Contours are plotted, and the 
density peak is marked with a circle. The colors of the 
peak are then compared with a BC06 track to obtain an 
age (shown in left panels of Figure [15]). 

4.1.1. NGC 520 

The interaction history of NGC 520 is the least well 
understood of our sample. Stanford & Balcells (1991) 
ran the most recent simulation of the system, comparing 
to infrared and optical images, as well the H I velocity 
field. Their best match to these observations are at a 
simulated age of ^ 300 Myr ago. The system is particu¬ 
larly complicated, however, because Stanford & Balcells 
find that their simulations were not fully consistent with 
two colliding disk galaxies, due to the presence of a third 
dwarf galaxy, UGC 957. They conclude that the south¬ 
ern tail formed from the interaction of two massive disk 
galaxies and that the northern tail is the result of a pas¬ 
sage by UGG 957. If the two tails formed via different 
interactions, we might expect differences in their cluster 
ages. Despite these complexities, we adopt a best guess 
of Ttaii ^ 300 Myr. 

The northern tail of NGC 520 contains only three 
young clusters, making a color density plot impractical 
here. However, the three clusters are tightly grouped 
around 100 Myr on the BC06 track; we therefore approx¬ 
imate their ages to be Tdus ^ 100 Myr. The majority of 
the south tail clusters appear to be older, and, based on 
their color distribution in Figure [TS] they are consistent 
with a period of cluster formation that began around 
Tdus ^ 260 Myr ago. 

The differences in the numbers and ages of the clus¬ 
ters detected in the north and south tails appear to be 
consistent with the scenario put forth by Stanford & Bal¬ 
cells. The passage of a dwarf galaxy would not necessar¬ 
ily result in a starburst as intense as one resulting from 
two equal-mass galaxies, because the dwarf galaxy would 
strip less gas to form clusters. We detect many more 
clusters in the southern tail, which Stanford & Balcells 
predict is the result of a disk-disk interaction, than the 
northern tail, which likely formed from the passage of the 
dwarf galaxy UGG 957. Stanford & Balcells speculate 
that the dwarf galaxy might have formed the northern 
tail as the disks were interacting, around the same time 
the southern tail was forming. However, the tight age 
grouping of the three clusters in NGC 520N seems to in¬ 
dicate that cluster formation took place earlier than in 
NGC 520S. 

4.1.2. NGC 2623 

Privon et al. (2013) simulated NGC 2623 and esti¬ 
mated that an age rtaii = 220 ± 30 Myr best reproduces 
spatial and kinematic distributions of H I. Privon et al. 
also successfully recreated the pie wedge region, showing 
that it is probably material from the northern tail that 
has fallen back through the main body and is currently 
on a southward trajectory. 


We measure clusters in NGC 2623N to have a peak 
density in color corresponding to an age of Tdus ^ 230 
Myr; this supports the scenario that cluster formation 
began at the same time that the tail formed. There are 
some clusters that appear to be older than the tail, how¬ 
ever, as we discussed in Section 3.3.1, we believe this is 
due to photometric scatter resulting from low S/N. NGC 
2623S possesses too few clusters to reliably measure Tdus, 
although it is clear from Figure [9] that most of the clus¬ 
ters are slightly bluer than the color associated with the 
tail age. It is possible that formation of clusters in the 
south tail was somewhat delayed relative to the tail itself, 
but uncertainties in tail and cluster ages are too large to 
definitively establish such a delay. 

Inspection of Figures [9] and [15] shows that clusters in 
the pie wedge have a very narrow range of colors, sug¬ 
gesting that they all formed at approximately the same 
time; we measure Tdus to be ^ 100 Myr. These clus¬ 
ters are clearly bluer and hence younger than clusters 
in the tails of NGC 2623. This result is consistent with 
the scenario proposed by Evans et al. (2008), where the 
pie wedge formed from debris from the inner region of 
the northern tail. They suggest that the debris fell back 
through the main body, inducing a burst of star forma¬ 
tion (see also Privon et al. 2013). While Privon et al. 
do not provide an age for the pie wedge, their Figure 
8 shows snapshots of the merger at different times. At 
85 Myr ago, the debris can be seen falling southward, 
past the main body, roughly consistent with our cluster 
formation timescale. 

4.1.3. NGC 3256 

Recent unpublished results from simulations place 
NGC 3256 at rtaii = 450 ± 50 Myr (G. Privon, private 
communication). The simulation is based on spatial and 
kinematic H I images of the system and was run using 
the same methodology as for NGC 2623. 

Clusters in the eastern tail have the largest spread in 
color of any system in our sample, with a weak density 
peak, Tdus ^ 260 Myr. This suggests a longer duration of 
cluster formation in this tail. These ages are consistent 
with previous results. Tranche et al. (2007) obtained 
spectroscopic ages for three clusters in the western tail, 
and found two clusters with ages of ^ 80 Myr, and a 
third with an estimated age of ^ 230 Myr. The latter 
fits well within the range that we find for clusters in the 
eastern tail. 

4.1.4. The Timescale of Cluster Formation 

In the three systems studied here, we have found that 
cluster formation typically begins with the tidal tail for¬ 
mation, although in a few cases a delay cannot be ruled 
out. This type of delay was also found by Bastian et 
al. (2005) in NGC 6872. We therefore conclude that, in 
general, Tdus < 'Ttaii- 

r < 10 Myr clusters have not formed in any of the 
tails studied here. We confirm this by examining HST 
Ha images of the two most recent mergers, NGC 2623 
and NGC 520, and find no Ha emission in any of their 
tails. We speculate that either the tails have used up 
their gas, or, because the gas is not collisionless, most of 
the remaining gas was largely dispersed during the inter¬ 
action. It is worth noting, however, that clusters with 
age r < 10 Myr have been reported in the tidal tails of 
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Figure 15. Color-color diagrams used to quantify when the peak cluster formation occurred in each system where possible. Right panels 
show densities of clusters in color-color space, with contours overplotted. The peak is marked with a blue circle. Left panels are similar to 
Figure [Q] with density peaks from the right panels included. 


the Tadpole galaxy and NGC 6872 (Tran et al. 2003; 
Bastian et al. 2005). In the case of NGC 6872, the time 
since tail formation is ^ 145 Myr (Horellou & Koribal- 
ski 2003), making NGC 6872 significantly younger than 
any interacting system in our sample (except for the pie 
wedge). It is possible that NGC 6872 has simply not 
yet dispersed its tidal tail gas, resulting in newly formed 
clusters. 

Based on our cluster age distributions, we suggest the 
following scenario: initially, the progenitor galaxies form 
clusters at a lower, approximately constant rate. Once 
the two galaxies begin their interaction, many of these 
disk clusters, along with other disk material, are stripped 
and form the tidal tails. The stripped gas typically be¬ 
gins forming clusters immediately, resulting in an in¬ 
creased rate of cluster formation. The process contin¬ 
ues in the tails at a lower rate until the gas reservoir is 
exhausted or the gas is dispersed. 

4.2. Comparing Cluster Ages Aeross Galaxies 

We suggest the following sequence of merger ages: 
NGC 2623 is the most recent of the three mergers, fol¬ 
lowed by NGC 520, with NGC 3256 being the oldest. 
This sequence is supported by both our estimated clus¬ 
ter ages and simulated ages of tidal tails (and thus ages of 
the mergers themselves). We note that the age sequence 
suggested by these simulations is different than what the 
Toomre Sequence suggests, as NGC 520 and NGC 2623 
are flipped (Toomre 1977). 


We find that clusters in the pie wedge have ages around 
R;ius ^ 100 Myr. The north tail of NGC 2623 contains 
clusters with ages at Tdus ^ 230 Myr, while the few 
south tail clusters are slightly younger, although we do 
not attempt to quantify their ages. NGC 520, the next 
youngest merger, exhibits a slight age difference among 
clusters in different tails. NGC 520N has clusters around 
R:ius ^ 100 Myr, while NGC 520S clusters has Tdus ^ 260 
Myr. The oldest merger, NGC 3256, contains clusters in 
the eastern tail with ages Tdus ^ 260 Myr, although we 
caution that the broad age distribution in NGC 3256E 
suggests that cluster formation likely began earlier than 
260 Myr ago. For a direct comparison, all the values of 
Ttaii and Tdus are listed in Table [H 

4.3. How are the Tidal Tail Age and Luminosity 
Funetion Related? 

We have measured the LFs for all tails where possible, 
and find tentative evidence that younger mergers exhibit 
shallower LFs, although our sample size is small. Mil 
found that a LF with —2.5 < a < —2 provided quali¬ 
tatively acceptable fits to their statistically determined 
cluster sample in tails where the measurement was pos¬ 
sible, although they were not able to directly fit the LF 
with a power law due to the poor quality of the data. 
Whitmore et al. (2014) measured the luminosity func¬ 
tion in a variety of spiral galaxies in search of various cor¬ 
relations between a and galactic environment. A weak 
negative correlation was found between star formation 
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rate (SFR) and a (i.e. galaxies with high SFR tend to 
have flatter LFs). For instance, the Antennae, a major 
merger with a relatively high SFR, had the second lowest 
a value in their sample, at a = —2.07 ± 0.03. 

No measurements of the cluster LF were made in tidal 
tails, however. While a tends to decrease with the tail 
age, all a values in our sample are within the range ob¬ 
served in a wide variety of other environments (Whitmore 
et al. 2014 and references therein). 

5. CONCLUSIONS 

We have used ACS/WFC observations from HST to 
directly observe star clusters in the tidal tails of three 
nearby mergers. We draw the following conclusions. 

1. Every tidal tail in our sample, as well as the pie 
wedge, contains a population of star clusters. We 
note that NGC 520N and NGC 3256E were re¬ 
ported to have no statistically significant cluster 
presence in Mil. 

2. We estimated the ages of clusters in each tidal tail 
and compared with estimated ages for the tails 
themselves from simulations. The tidal tails in 
NGG 2623 contain a population of clusters which 
appear to have formed during the formation of the 
tails. Tails of NGG 520 and NGG 3256 possess sev¬ 
eral clusters that appear to predate the merger, in 
addition to a population of clusters that formed at, 
or soon after the formation of the tails. 

3. In every tidal tail in our sample, cluster formation 
lasted for several tens of millions of years, but no 
very recent formation (in the last few million years) 
has occurred in any of the tails. 

4. Simulations place the formation of tails in the fol¬ 
lowing sequence, from youngest to oldest: the pie 
wedge, NGG 2623, NGG 520, NGG 3256. Ages 
obtained from clusters generally agree with this se¬ 
quence. 

5. The diffuse light in the tails of NGG 520 and NGG 
2623 exhibits a gradient in color across the tail (as 
opposed to along it), which is likely indicative of 
systematic patterns in the ages of the diffuse stel¬ 
lar population. To our knowledge, such color gradi¬ 
ents in the diffuse stellar light in tidal tails of these 
galaxies has not been previously reported. The gra¬ 
dient is loosely traced by the spatial distribution of 
young and old star clusters within the tail. We in¬ 
terpret this gradient as a superposition of broadly 
distributed older stellar material and younger stars 
and clusters that formed along the center of the 
tail, where gas is densest. 

6. The LEs of clusters in our tidal tail sample are sim¬ 
ilar to those found in a variety of galactic environ¬ 
ments. We tentatively find that as the merger age 


increases, the luminosity function tends to become 
steeper. 
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